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Abstract. We show in the smooth category that the heat trace asymptotics 
and the heat content asymptotics can be made to grow arbitrarily rapidly. In 
the real analytic context, however, this is not true and we establish universal 
bounds on their growth. 
MSG 2002: 35K20,35P99,58J25,58J50 



1. Introduction 



, 1.1. Heat trace asymptotics. Let {M,g) be a compact Riemannian manifold of 

' dimension m with smooth (possibly empty) boundary dM . Let dvolm and dvol,„_i 

be the Riemannian volume elements on AI and on dM, respectively. Let Ag be 
the scalar Laplacian. Let v be the inward unit normal on the boundary; we extend 
ly by parallel translation to a vector field defined on a collared neighborhood of 
the boundary so Vj/J^ = 0; this means that the integral curves of v are unit speed 
geodesies perpendicular to dM. Let 

^ . B^c/) := (j>\dM and B^(f) i^(/)|aA/ 

be the Dirichlet and Neumann boundary operators, respectively. Impose boundary 
(y-^ I conditions B = B^ or B — B^ . Let u : M x (0, oo) — ^ R be the unique solution of 

{dt + Ag)u(.x,t) — (heat equation), 

limt_>.o u(-; t) 4>i{') in (initial condition), 
Bu{-, i) = for i > (boundary condition), 

where 4>i is real-valued and smooth on M . Then u{x, t) represents the temperature 
^ . at a; € M at time t > ii M has initial temperature distribution where the 

boundary condition B is imposed on u for t > 0. The solution is formally given by 

where Ag g is the associated realization of the Laplacian. The operator e^^^"-'^ is 
a smoothing operator of trace class and, as < J, 0, there is a complete asymptotic 
series of the form [29, 30, 43, 44, 45, 46, 47] 

OO 

Tri.{e-*^-n ^ (4^t)-™/2^a„(M,5,S)i"/2. 

If M is a closed manifold, the boundary condition plays no role and we shall denote 
these coefficients by a„(-M, g). They vanish if n is odd in this instance. 

The asymptotic coefficients {ai, a2, • ■ ■ } are locally computable invariants of M 
and of dM as we shall see presently in Section 2. In mathematical physics, they 
occur for example in the calculation of Casimir forces [5, 18, 33] or in the study of 
the partition function of quantum mechanical systems [6, 8, 33]. They are known 
in the category of manifolds with boundary for n < 5 [19, 32], and in the category 
of closed manifolds for n < 8 [1, 4]. These coefficients play a crucial role in the 
study of isospectral questions. Related invariants for more general operators of 
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Laplace type also play a crucial role in the local index theorem. See, for example, 
the discussion and references in [2, 3, 20, 21, 23, 26, 27, 28, 37, 38, 39]. They have 
also been studied with nonlocal boundary conditions [34]. We also refer to [24] 
where the heat trace itself is studied and not just the asymptotic coefficients. For 
the study of the asymptotic behaviour of the eigenvalues of ^g,B we refer to [41] 
and the references therein. The field is vast and it is only possible to cite a few 
references. 

1.2. Planar domains. In the case of a planar domain f2, the heat trace asymp- 
totics (with Dirichlet boundary conditions) have been computed for n < 13 by 

Berry and Howls [17]. Berry and Howls computed a„ for n < 31 in the case of a 
disc [17], and were led to conjecture that for planar domains ^l and for n — )■ oo 

an{^) = aV{n - /3 + l)T{n/2)-^£{Q.f-''{l + o(l)), (l.a) 

where a and (5 are dimensionless quantities and where is the length of the 
shortest accessible periodic geodesic in VI. In particular, for a disk of radius R and 
shortest accessible periodic geodesic 4i?, they further conjectured that Equation 
(l.a) holds with a = (8-\/27r)~^ and /? = |. While the latter conjecture remains 
open to date, it is instructive to see that Equation (l.a) can not hold in general. 
The following counter examples were given in [7] . 

Example 1.1. Let < e < ^, and let 

Pe = {{xi,X2) e K2 : < 1, 1.T2I < 1 - e}, 

Qe = {{xi,X2) e : < l,xi < 1 - £,2:2 < 1 - £}. 

We smooth out the corners of dPg at 0:2 = ±(1 — e) and of OQ^ at xi = 1 — e, 
X2 = I — s isometrically to obtain two convex domains and with smooth 
boundary and with o„(Pe) a„(Qe) and £(Pe) = 4(1 - e), i{Qs) = 2(2 - e). This 
then contradicts Equation (l.a). 

Example 1.2. Let < £ < 1, < p < 1 - £, and let 

D.^ := {{xi,X2) eE.^ ■.e<\x\< 1}, 

nP := {(.Ti,a;2) G . |^| < _ > 

We then have that an{^e) = an{^l) and = 2(1 - s), 1{VLP) = 2(1 - £ - p) 

which once again contradicts Equation (l.a). 

It remains an open problem to construct a pair of iso - a„ real analytic simply 

connected planar domains which have different shortest periodic geodesies. It has 
been conjectured that Equation (l.a) also holds for balls in R™ where j3 depends 
upon d only [31]. 

1.3. The heat trace asymptotics in the real analytic category. The calculus 

of Sceley [43, 44, 45, 46, 47] and Greiner [29, 30] shows that a„ is given by a local 
formula; the following result will then follow from the analysis of Section 2: 

Theorem 1.1. Let B be either Dirichlet or Neumann boundary conditions. There 
exist universal constants Kn,m so that if (M, g) is any compact real analytic manifold 
of dimension m, then there exists a positive constant C = C{M, g) such that 

\an{M,g,B)\ < C" • volm(M, g) for any n . 

We note some similarity between the formulae of Equation (l.a) and Theo- 
rem 1.1. The geometric data of {M,g) appear in C", whereas the prefactor is of a 
combinatorial nature and depends on m and n only. We can choose the constant 
to rescale appropriately under homotheties, i.e. so that C{M,c^g) = c~^C{M,g). 
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We restrict momentarily to the context of closed manifolds, i.e. compact man- 
ifolds with empty boundary. We adopt the Einstein convention and sum over re- 
peated indices. We say that D is an operator of Laplace type, if in any local system 
of coordinates we may express D in the form: 

D = - (g^'d^^d,^ + A'^d,, + B) . (l.b) 

Let an{x,D) be the local heat trace invariant of such an operator. We shall pri- 
marily interested in the case n even so we shall set n = 2n in what follows. If / is 
any smooth function on M, then 

Tr^=(/e-*'') ~ (47rt)-™/2 V i" / a^nix, D)f{x) dvol„ . (l.c) 

The following result shows that the factorial growth conjectured by Berry and Howls 
for planar domains pertains in this setting as well as regards the local heat trace 
invariants on closed manifolds. 

Theorem 1.2. Let {M,g) be a closed real analytic Riemannian manifold of dimen- 
sion m > 2. 

(1) Let D be a scalar real analytic operator of Laplace type on M. Then there 
exists a constant Ci = Ci{M,g,D) so that 

\a2n{x, D)\ < C" • n! for any n > 1 . 

(2) Let P be a point of M . Suppose there exists a real analytic function f on M 
such that df{P) ^ 0. Then there exists a constant C'2 = C2{P, M,g, /) > 
and there exists a real analytic function h on M so that the conformally 
equivalent metric gh ■= e?'^g satisfies 

I a2n (P, Ag J I > C2" • n! for any n > 3 . 

Remark 1.1. Assertion (1) can be integrated to yield an upper bound on the heat 
trace asymptotics a2n{D)- However, Assertion (2) is only valid at a single point of 
M. Since it in fact arises from considering a divergence term in the local expansion, 
we do not obtain a corresponding estimate for a2n{D). 

1.4. The heat trace asymptotics in the smooth category. The situation in 
the smooth non real analytic setting is very different. Fix a background reference 
Riemannian metric h and let V'* be the associated Lcvi-Civita connection which 
we use to covariantly differentiate tensors of all types. If T is a tensor field on M, 
we define the C'' norm of T by setting: 

Changing h replaces ||T||fc by an equivalent norm; we therefore suppress the depen- 
dence upon h. But as we will be changing the metric when considering the heat 
trace asymptotics subsequently, it is useful to have fixed h once and for all so the 
associated norms do not change. Theorem 1.1 fails in the smooth context as we 
have: 

Theorem 1.3. Let k > 3 be given, let constants Cn > for n> k be given, and let 
e > be given. Let {M, g) be a smooth compact Riemannian manifold of dimension 
m>2 without boundary and let g^ be the usual Euclidean metric on R"*+^. 

(1) There exists a function f e C°°{M) with ||/||/s_i < e so that if gi := e^^ g 
is the conformally related metric, then 

\a2n{.M , g-i)\ > Cn for any n>k. 
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(2) Suppose that g = &*ge where 8 is an immersion of M into ffi.™"*"^. There 
exists an immersion ©i with ||0 — Gi||fe-i < e so that if gi := @*ge, then 

\a2n{M,gi)\ > Cn for any n>k. 

1.5. Heat content asymptotics. There are analogous results for the heat content 

asymptotics. Let be the initial temperature of the manifold and let 02 be the 
specific heat of the manifold. We suppose throughout that and (1)2 are smooth. 
The total heat energy content of the manifold is then given by: 

l3{(pl,(j)2,Ag,B){t):= / u{x,t)(p2{x)dYo\m ■ 

Jm 

As 1 10, there is a complete asymptotic expansion of the form 

p{ci>i,4>2,Ag,B)it) ~ y2—r ^"gh-hdvoU 

oo 

+ j2t^'+'y'pr{^„^2,Ag,B). 

The coefficients involving integrals over M arise from the heat redistribution on 
the interior of the manifold and are well understood. The additional boundary 
terms (3f^ are the focus of our inquiry. They, like the heat trace asymptotics, 
are given by local formulae and have been studied extensively (see, for example 
[12, 13, 14, 15, 16, 22, 35, 36, 40, 42] and the references contained therein). 

Inspired by the work of Howls and Berry [31], Travenec and Samaj [48] inves- 
tigated the asymptotic behaviour of the coefficients as £ — >■ cxd in flat space in 
the special case that = <;A2 = 1 with Dirichlet boundary conditions. The inte- 
rior invariants then play no role for n > 1 and one has, adopting the notational 
conventions of this paper, that 

oo 

^(1, 1, Ag, B-){t) ~ vol™(M, g) + Y, t^'+'^/'^f ^(1, 1, A„ B") . 

1=0 

After interpreting the results of [48] in our notation, they found that if M is a ball 
in of radius r with m even, then as £ — )• oo one has: 

= 47r(™-3)/2p(^/2)-i(£ + l)-ir(£/2)r'"-^-i(l + o(l)) . (l.d) 

The structure of Equation (l.d) is similar to that of Equation (l.a). There is a 
combinatorial coefficient in m and f , while the shortest periodic geodesic appears 
to a suitable power. However, for m odd Travenec and Samaj obtained polynomial 
dependence rather than factorial dependence of in ^ [48] . Furthermore the two 
examples in Section 1.2 above provide iso-/3^ pairs of smooth planar domains with 
different shortest periodic geodesic lengths. Hence the structure of the asymptotic 
behaviour of the ^^'s in flat space remains unclear in general. 

For ^ even, the boundary term involves a fractional power of t and there is no cor- 
responding interior term. This simplifies the control of these terms. Consequently, 
we shall usually set ^ = 2^ in what follows. 

1.6. The heat content asymptotics in the real analytic setting. As noted 
above, results of [48] showed that the heat content asymptotics on the ball in M™ 

for m even exhibit growth rates similar to that given in Theorem 1.2 for the local 
heat trace asymptotics. We generalize Theorem 1.2 (2) to this setting to derive an 
estimate using conformal variations which shows that the metric on the boundary 
does not play a central role in the analysis: 
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Theorem 1.4. Let m>2. 

(1) Let {N,g]\r) be a closed Riemannian manifold of dimension m — 1. Let 

M := [0,2tt] X A'^. There exists a real analytic function h{x) on [0, 27r], 
which depends on the choice of {N^qn), so that the conformally adjusted 
metric qm ■= e^^{dx'^ + qn} satisfies: 

|/3ff (l,l,Ag„,B-)| >Z!-voU_i(Ar,5^) for any £>3. 

(2) Let Qf, he the standard Euclidean metric on the unit disk -D™ in . There 

exists a radial real analytic function h on , which depends on m, so that 
the conformally adjusted product metric gjyj :— e^^g^. satisfies: 

|/3ff (1, 1, Ag,, ,B-)\> i\ ■ vol„_i (iV, gN) for any l>'i. 

We have estimates for the heat content asymptotics in this setting which are 
similar to those given in Theorem 1.1: 

Theorem 1.5. There exist universal constants Kn,m and Ke,m such that if {M,g) 

is a compact real analytic Riemannian manifold of dimension m and if {(pi, 4)2) are 
real analytic, then there exists a positive constant C = C{M,g,(l)i,(j)2,B) such that 



f 



1 • A^(/)2 dvol„ 

M 



< K„,mC" -vol^ (M,c/), 



Remark 1.2. Again, the constant C can be chosen so that 

CiM,c-'g) = c-^C{M,g). 

1.7. The heat content asymptotics in the smooth setting. Theorem 1.5 fails 

in the smooth setting as wc have: 

Theorem 1.6. Let k > 3 be given, let constants Cf > for £ > k be given, and 
let e > be given. Let B = B+ or B = B~ . Let (M, g) be a smooth compact 

Riemannian m,a,nifold of dimension m > I with non-trivial boundary. Let cpi be a 
smooth initial temperature and let <p2 be a smooth specific heat with B(j)2 ^ 0. There 
exists $1 with \ — $i||2fc-i < £ such that: 

f4f{^u<t>i,^g,B)=Ci for any l>k. 

The heat content asymptotics were originally studied for Dirichlet boundary 
conditions and for (pi = (p2 = 1 [9, 10, 11]. We have the following theorem in this 

setting: 

Theorem 1.7. Let k > 3 be given, let constants Ci > for I > k be given, 
and let e > be given. Let {M,g) be a smooth compact manifold Riemannian 
manifold of dimension m > 2 with non-trivial boundary. There exists a metric gi 
so Wq — 9i\ \2k-1 < £ such that 

/3ff (l,l,Ag,,S-) = Q- for any £>k. 

1.8. Bochner formalism for operators of Laplace type. The results given 
above in Theorem 1.3, in Theorem 1.6, and in Theorem 1.7 rely upon a leading 

term analysis of the heat trace asymptotics and of the heat content asymptotics. 
It is one of the paradoxes of this subject that to apply the functorial method, one 
must work with very general operators even if one is only interested in the scalar 
Laplacian, as is the case in this paper. Wc only consider the context of scalar 
operators. There is a corresponding notion for systems, i.e. operators which act 
on the space of smooth sections to some vector bundle. It is possible to express an 
operator D of Laplace type as given in Equation (l.b) invariantly using a Bochner 
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formalism [27]. There exists a unique connection V and a unique smooth function 
E so that 

where we use ';' to denote the components of multiple covariant differentiation 
with respect to V and with respect to the Levi-Civita connection. Let T^v^ be 
the ChristofFcl symbols of the Levi-Civita connection and let ui be the connection 

1- form of V. We then have 

E = B-g''''{d^^uj^+uJuOJv-^w^uv'")- 

1.9. Leading term analysis. Theorem 1.8 below will play a central role in our 
analysis, and was established in [20, 25, 26]. We also refer to related work in the 

2- dimensional setting [39]. It has been used by Brooks, Perry, Yang [21] and by 
Chang and Yang [23] to show families of isospectral metrics within a conformal 
class are compact modulo gauge equivalence in dimension 3. Let t be the scalar 
curvature of g, let p be the Ricci tensor of and let i7 be the curvature of the 
connection V defined by an operator of Laplace type. 

Theorem 1.8. Let D he an operator of Laplace type on a closed Riemannian 

manifold (M,g) and let n> 3. 

(1) The local heat trace asymptotics satisfy: 

a2n{P,A,) = ^-L-^{-nA^-^T-{4n + 2)A^-'E} 

+ lower order derivative terms. 

(2) The global heat trace asymptotics satisfy: 

a.n{D) = i^^/{(n2-n-l)|V^-V+2lV^-V 
2 (2n+ Ij! Jm 

+ 4(2n + l)(n - 1)V("-2V • V^^-'^^E + 2(2fi + 1)] V^"-^)^]^ 
+ 4(2n - l)(2n + 1)] V""^£;]^ + lower order terms } dvol^ . 

In this paper, we will establish a corresponding leading term analysis for the heat 
content asymptotics. We shall always assume £ is even; thus the lack of symmetry 
in the way we have written the interior contributions plays no role. Let V be the 
connection defined by D as discussed in Section 1.8. Let D* be the formal adjoint 
of D\ the associated connection V* defined by D* is then the connection dual to 
V defined by the relation 

V01 • + f^l • V*(^2 = d{(j)i ■ (j)2) ■ 

Let 

(/.f) := V5</.i|9M and 4'^ := (V:)V2|9m 
be the normal covariant derivatives of order i. By using the inward geodesic flow, we 
can always choose coordinates {y, r) near the boundary so that dr = v, consequently 

ct>^^^ = di<k\dM if D = Ag. 

Let 5 be a smooth function on the boundary. The Robin boundary operator in this 
more general setting is defined by the identity: 

DM ■ 

Let Pram '■ — R 

am7na\7n . . .m t>e tlie covariaiit derivative of pmm 

restricted to dM. 

Define recursively for i even by setting: 

S2 = -2^-i/2| and S^ = j^S^_2 if ^>4. 
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Theorem 1.9. Let £ > 6 be even. Modulo lower order terms we have: 



.IrtM 



JdM 



+0 • ^('V^'^^^^-*' + - 2)E(<p,<P2p'L?^ + ... } dvoWi. 



. flM 



JdM 



+(2 - e)E,4>^^U^2^E('~^) - E,Si<pt'U2 + Mt'^) 



Pmm 



1.10. Outline of the paper. In Section 2 wc will prove Theorem 1.1 and The- 
orem 1.5. In Section 3, we use Theorem 1.9 to establish Theorem 1.6 and Theo- 
rem 1.7. In Section 4, we use Theorem 1.8 to demonstrate Theorem 1.3. Theo- 
rem 1.9 is new and is proved in Section 5 by extending functorial methods employed 
in [12, 13]. In Section 6, we establish Theorem 1.2. We conclude the paper in Sec- 
tion 7 by demonstrating Theorem 1.4. 



\a\:=ai-\ h ■= {dxiT' ■ ■ ■ {dx^T"' , 9i]/a — d'^9t] for |a| > 0. 



In any local system of coordinates, the Riemannian volume form on M is given by: 



Let g^^ be the inverse matrix; this gives the components of the dual metric on the 
cotangent bundle. Since the heat trace and heat content asymptotics are given by 
suitable local formulae, Theorem 1.1 and Theorem 1.5 will follow from the following 
result: 

Theorem 2.1. Let £n he a local interior invariant which is homogeneous of degree n 
in the jets of the metric and a finite (possibly empty) collection {0i, ...} of additional 
smooth functions. Let Tn-i be a local boundary invariant which is homogeneous of 
degree n—1 in the jets of the metric and a finite (possibly empty) collection {4>i, ...} 
of additional smooth functions. Let (M, g) be a compact real analytic manifold of 
dimension m with real analytic (possibly empty) boundary dM so that the metric 
g is real analytic and so that the collection {0i, ...} is real analytic. There exists a 
constant C = C(M, > (which is independent of the choice of and of 

Tn) and there exist constants n{£n) > and K{Tn-i) > (which are independent 
of the choice of {M,g,(j)i, ...)) so that 



2. Local invariants in the real analytic setting 



Let a := (ai, . . . , 



Olm) be a non-trivial multi-index. We define: 




(x, 5, </>!,. ..)dvolm < K{£n)C ■^o\m{M,g) 



JM 




The constant C{M, g, (pi, ...) may be chosen so that 

C{M,c^g,<Pi,...)=c-"C{M,g,^i,...). 
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Proof. Suppose first that the boundary of M is empty. For each point P of M, there 
exists e{P) > so the exponential map defines a real analytic geodesic coordinate 
ball of radius e{P) about P. Let /C be a compact neighborhood of the identity in 
the space of all symmetric m x m matrices. Since gij = 6ij at the center of such 
a geodesic coordinate ball, by shrinking e{P) if necessary, we may assume that the 
matrix (gij) belongs to JC for any point of the coordinate ball of radius e{P). Since 
we are working in the real analytic category and since {gij, <pi, ■■■} are real analytic 
near P there exists a C = C{P,M,g,(j)i, ...) so that again by shrinking e{P) if 
necessary we have that 

\dZg^j\<C\'^^a\\ and MS</.^| < ^"1 |a|! on S,(p)(P) (2.a) 

for any miilti-indcx a. Wc cover M by a finite mimbcr of such coordinate balls 
about points (Pi,...) and set C(M, 5, ...) = maxj, C(P^, M, t;, 0i, ...). Since £ is 
a local invariant, we may expand: 

e{x,g) =Y.e.^^^{gij{x)m^gi,j,)...{d:^gi^j^) ■ {4^M...{4'<i>k,) (2.b) 

where in this sum we have the relations: 

|ai| + ...\aa\ + + ... + \Pb\ =n, < \ai\, < |q!o| • 

Since e - ^ is continuous on the compact neighborhood /C of the identity 6, we may 
bound 

\^S.0(9tj{x))\ < E^ uniformly on /C . 

Combining the estimates of Equation (2. a) with the estimates given above and 
summing over (a, /?) in Equation (2.b) yields an estimate of the desired form after 
integration. Since f„ is homogeneous of degree n, it follows that 

£n{x,c'^g,^i,:.) = c~'^£n{x,g,(j)i,...) . 

The desired rescaling behaviour of the constant C{M,g,(pi, ...) now follows. 

If the boundary of M is non-empty, we must also choose suitable coordinate 
charts near dM. If Q e dM, we consider the geodesic ball B^^{Q) of radius e in 
dM about Q relative to the restriction of the metric to the boundary and wc shall 
let B^^,{Q) := [0, t) x B^^^^{Q) for some i > be defined using the inward geodesic 
flow so that the curves r (r, Q) are unit speed geodesies perpendicular to the 
boundary. Again, by shrinking e and l, we may achieve the estimates of Equation 
(2. a) uniformly on B^^, {Q). We cover M by a finite number of coordinate charts 
Bs{P) for P e int(Af) and B.^^iQ) for Q G dM . The desired estimate for now 
follows. To study the invariant Tn-i, we cover dM by a finite number of coordinate 
charts B^^^iQ) for Q S dM and argue as above. □ 

3. Leading Terms in the Heat Content Asymptotics 
We shall omit the proof of the following result as it is well known. 
Lemma 3.1. 

(1) Let k > 1 be given, let constants 7^ > for £ > k be given, and let e > 
be given. Let {M,g) be a smooth Riemannian manifold with non-empty 
boundary dM. There exists a smooth function ^ on M so that \ < e 
and so that 

ci)W = '4'iy)^e for £> k. 

(2) Let k>l be given, let C > be given, and let e> be given. There exists 
a smooth function f on M := [0, 1] with ||/||fc_i < e and Jj^ \d^f\'^dx > C. 
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The proof of Theorem 1.6 and of Theorem 1.7. Let fc > 3 be given, let constants 
> for £ > fc be given, and let e > be given. Let (M, g) be a smooth compact 
Riemannian manifold of dimension m>2 with non-trivial boundary. We first take 
B = B~ to consider Dirichlet boundary conditions. Let 4>i be a smooth initial 
temperature and let 02 be a smooth specific heat with B~02 7^ 0. Since (/>2 does 
not vanish identically on the boundary, there exists a smooth function ij) on dM so 

%l}4>2 dvol^-i = 1 . 

dM 

Let {71, ...} be a sequence of constants, to be determined presently. For v > k, let 



r 



2j 



j=k 

Since /32/ is given by a local formula of degree 2£, only the constants 71, 7^ play 
a role in the computation of (3^^^ , i.e. 

/3ff + </>2,Ag,S) = /3ff ($^- +</)!, (/.2,Ag,S) if fi>£. 

We take ^k-i = 0. Since 7^ 0, we can recursively choose the constants 7^, 
and hence the functions for £ > A; so 

^21 ■ Te = C-e - + cPuh, ^g,B) for £>k 

and apply Theorem 1.9 to see: 

Plfi^2l + <t>iA2,A„B-) = Ce. 

We complete the proof of Theorem 1.6 (1) by using Lemma 3.1 to choose $ with 
||*||2fc-i < £ such that 

if j < 2k or if j is odd 
7^- if j = 2£ for i>k 

To prove Assertion (2) of Theorem 1.6, we use Assertion (2) of Theorem 1.9 and 

examine the term — S2£(f>i^^ ^^2^'' ^ to prove Theorem 1.7, we apply Assertion (1) of 

Theorem 1.9 and examine the term i(2£ — 2)^2l4'i4'2Pmm As apart from these 
minor changes the proof is exactly the same as that given above, we shall omit 
details in the intersts of brevity. □ 

4. Leading terms in the heat trace asymptotics 

4.1. Proof of Theorem 1.3 (1). We set = and = in Theorem 1.8 to 
study the Laplacian and see thereby that there exists a non-zero constant d„ so: 



$0) 



• ' ,2 



a2n{Ag) = dn [ {(n2-n-l)|V"-V|2 + 2|V 

+ Qn,m(-R,VJ?,...,V"-3i?)}dvol™ . 

Let £ > be given. We restrict to a single geodesic ball B of radius 35 for some 

6 > about a point P. Let 6* be a plateau function so that 6* = 1 for |x| < 5 
and = for \x\ > 26. We shall define the functions fk, /fe+i, ••• recursively and 
consider the conformal deformation: 

g9(x)(2/fc(xi)+...+2/^(xi))^^ 

Let fc > 3. Choose < 5^ for fc < so that < 5^ for fc < /x implies: 

Constraint 4.1. 

(1) /oo := hm^^ool/fc H 1- fp,} converges in the topology for any £. 

(2) (7oo := hm^_>.oo converges in the topology for any £. 
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(3) ll/IU-i <£. 

(4) ||5;x-5/i+ilU < 2"''£ /or ant/ /X. 

A-priori, one must consider jets of degree 2n in computing 02,-1 (Ag) (and in fact 
this is the case when considering the local heat asymptotic coefficients of Equa- 
tion (l.c)). However, by Theorem 1.8, only the jets of the metric to degree n play 
a role in the computation of the integrated invariants, a2n- 

Constraint 4.2. Choose < (5^ < 5^ for /c < /U so ||/^||^-i < 5^ for < /x implies: 

(1) |a2fx(Ag^_ J - a2n(AgJ| < 2"" for 3 < < n < /x- 

(2) |a2n(AgJ| - 1 < M^gJl for 3 < < n. 

The polynomial Qn,m{-) involves lower order derivatives of the metric. 

Constraint 4.3. Choose < 5^ < 5^ for A; < /x so that ||/^||^-i < for k < fj, 
implies there are constants = Cj^{fk,--- , f^-i) depending only on the choices 

made previously so 

MAgJ > Kl / {|2V''-S.r + (M'-/^-l)|V"-Vl'}dvol„-C; 
> \d,\ [ {\2W-\f}dvolm-Cl. 

JBs 

On Bs, the plateau function 6 is identically 1 and we have: 

From this it follows that 

V"^^T = (m — l)d^_^ffi + lower order terms . 
Since gij is in a compact neighborhood of 6ij, we may estimate: 

l|V"-VgJp(P) > = \d^Jf,\^ + lower order terms . (4.a) 

Constraint 4.4. Choose < (5^ < (5j^ for /c < /x where (5^ = 5'^{fk, ■ ■ ■ , 
depends on the choices made previously so that ||/,i||,i-i < 6'^ for k < ^ implies 
there are constants = C^(//s, . . . , depending only on the choices made 
previously so 

JBri JBri 

(1 (1 

Theorem 1.1 (1) now follows from Lemma 3.1 (2). We can choose recursively 
subject to the constraints given above so that is arbitrarily small and so 

that 1 9^^/^ I ^ dvol^ is arbitrarily large. □ 

4.2. The proof of Theorem 1.1 (2). Let {M,g) be a hypersurface in M™. We 

fix P G M. After applying a rigid body motion, we may assume that P = 
and that the normal to M at P is given by Cm+i '■= (0, . . . , 0, 1). Thus we may 
write M as a graph over the ball in in the form x — )• {x,fo{x)) where 
/o(P) = and dfo{P) = 0. Let 9 he a. plateau function which is 1 for \x\ < 6 
and for I x I > (5. We shall consider the perturbed hypersurface defined near P 
by x ^ (xjoix) + eix)ifkix) + ...)) where /^(P) = and d/^(P) = 0. This 
hypersurface agrees with the original hypersurface away from P. We shall need 
to establish an analogue of Equation (4. a). The remainder of the analysis will be 
similar to that performed in the proof of Theorem 1.1 (1), and will therefore be 
omitted. 
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Suppose we have a hypersurface in the form \E'(a;) := {x,F{x)) where F{0) = 
and dF{Q) = 0. Let Fi := dxiF, Fij := dx^dx^F, and so forth. We compute: 

**(^ccj = ei + FiCm+l, 

^jki = ^{FjkFi + FjiFk + FjkFi + FkiFj — FjiFk — F^iFj} = Fji-Fi, 

rjfc' = g''"FjkFn, 

Rijk = g'^"{FjkFin - FikFjn} + lower order terms, 

where the lower order terms are either 4*^ order in the 1-jets or linear in the 
2-jets and quadratic in the 1-jets. We suppose F = -F)j_i + where we set 

//n = £m cos{af^x'^) cos(6^a;2). 

T = 4e^a^6^{cos^(a^a;^) cos^(6^a;^) — sin^(a^a;^) sin^(6jua;^)} + . . . , 
IV'-Vp = As^aXl cos2(a^.xi) cos2(6^.t1) - sm^{a^x^) sm'{b^x^)\^ + . . . , 

where wo have omitted lower order terms cither involving or not multiplied by 
the appropriate power of a^6J^. To simplify matters, we suppose 6 = n and that 
and 6^ are non-zero integers. We use the fact that we are dealing with periodic 
functions to compute: 

/•TT />7r 

/ / I cos^{af^x^) cos^(6^a;^) — sm^{af^x^) sm^{hf_ix'^)\'^dx'^dx^ 

Jx^ =—7r Jx^ =—7r 

= «;'^;' T"" /^"^ \cos\x^)cos\x^)-sm''{x^)sin\x^)\^dx^dx^ 

J X^= — a^7T J x'^ — — bfj_7T 

— ^j^I^^^^^p^M / / I cos^(a;^) cos^(x^) — sin^(x^) sm^{x'^)\^dx^dx^ 

Jx^ — — TT Jx^ — — 7r 

We shall take = a'^, take large, and take S/j, appropriately small to complete 
the proof. □ 

5. Leading terms in the heat content asymptotics 

This section is devoted to the proof of Theorem 1.9. Let D be an operator of 
Laplace type on a compact smooth Riemannian manifold (Af, g) with non-empty 
boundary. We adopt the notation estabhshed in Section 1.8 and in Section 1.9. We 
shall always take S to be real in defining the Robin boundary operator. One then 
has the symmetry 

^{4>r,4,2,D,B){t) = p{4>2Ai,D\ B) it) . (5.a) 

If i is even, the lack of symmetry in the way we expressed the interior terms plays 
no role and thus Equation (5. a) yields: 

/3ff (<^i,<^2,D,i3) = /3ff (</.2,<Ai,£'*,B). (5.b) 

Let indices {a, b} range from 1 to m — 1 and index the tangential coordinates 
{y^, . . . ,7/™"^) in an adapted coordinate system such that dr is the inward unit 
geodesic normal. We then have 

ds'^ = gab{y, r)dy°' o dy'' + dr o dr . 

We define the second fundamental form by setting: 

Lab ■■= 9{^dy^dy^,dr) = -^drQab ■ 

Results of [12, 13] yield the following formulae which will form the starting point 
for our analysis: 
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Lemma 5.1. Adopt the notation established above. Then 

(1) /^f*^ 02, = -^/3^0i02dvoU_i. 

(2) Pl^^{<t>u4>2,D,B+s)= 0. 

(3) /3f^(<^i,<^2,i?,e-) =-^/g^{| (^'^2 +<^1<^2£; 

+ {j2^aa^bb - \LabLab - \pmm) ^l<p2} dvolm-i. 

(4) l3l^{cj,,,cl,2,D,B+)^ +5</.2)dvol™_i. 

We begin the proof of Theorem 1.9 by expressing Pf'^, modulo lower order 
terms, in terms of certain invariants involving maximal derivatives with unknown 
but universal coefficients; the symmetry of Equation (5.b) plays a crucial role in 
our analysis. Standard arguments (see [12]) show the coefficients in the following 
expressions are independent of the underlying dimension of the manifold: 

pr{4>uh,D,B-)= I |c,-,(</.fV2+</'i4'^)+c,-2(<Ar^v^'^+0i'vr'^) 

JdM 

+rj 4>i<j)'2pm~^^ + ...| dvolm_i, 

pr{^^,cp„D,B+)= [ {c+,(<^fV2+<^i<^f)+c+2(<^r'v^'^+4'vr'^) 

JdM 

+d+^S{<j)i<j>i''> + <I>['U2)E^^~^^ + d+,S<j>i<j>2E^'-^'> + r+<j)ihplL^^ 
+...| dvolm_i. 

We will determine all the coefficients except rf^g in what follows. Recall that 
S2 = -27r"^/^^ and Eg = ■ 

Lemma 5.2. Let £>4 be even. Let B = B~ or B = Bg. 

(1) Let D he self-adjoint with respect to the boundary conditions defined by B. 
IfB^i=0, then pf^{^i,^2,D,B) = ^A-2(<^f^ + ^^,-^2, AB)- 

(2) = Ee, = 0, = 0, and c^j = 

(3) = e^3 = 0, e^^ = 0, = ~'^i> ^'"''^ '''t = 0- 

(4) rf+i = d+2 = -Ee. 
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Proof. We follow [12] to derive Assertion (1) as follows. Let {A^, 4>fj,} be a complete 
spectral resolution of Db- Here {^^} is a complete orthonormal basis for L'^{M) of 
smooth functions with Dcf)^ = Xfi(pix and Bcf)^ = 0. Let 



JM 



JM 

be the associated Fourier coefficients. Then 

•DO 

If B(j)i = 0, then 

7^(D0i) = / D(bi- dvoU = / 01 • Dc^^ dvoU = Kll^i^i) ■ 

JM JM 

Consequently we have that: 
/3(Z)<^i,02,^,S)(i) 

^LiL / D"+Vi-<^2dvoU + 5^#+i)/2^f^(i?</.i,</.2,AB) 

n JM ,,_n 



n=0 "' -^^^ fc=0 

= E '^-*'''7^'(i?<^i)7^'(</'2) = E Ve-*^''7^'(</'i)7^'(</'2) 

/.i— 1 /J.— 1 

= -|Ee"*'''<('^i)^?('^2) = -|/3(<^i,<^2,AS)(t) 

~ I I?^0i-<^2dvoU-E^i(^-i)/^/3r(,/.i,</.„AB). 

The asymptotics defined by the interior integrals are the same. We note that 

(2) 

—D(j)i = + E(j>i. We set fc = £ — 2 and equate the asymptotics defined by the 
boundary integrals to establish Assertion (1). 

l{£ = 2, then the relations of Assertion (2) would follow from Lemma 5.1 modulo 
the caveat that we have but a single term c^2'P\^'t^2^ rather than 2 distinct terms 
in that setting. This will let us apply the recursion relation of Assertion (1) even 
if ^ = 4. Let 4'i\dM ~ 4>i\dM = 0. We set E = Q and consider cfj^(f)^p(j)2 and 
c^2</'i^~^V2^''- These terms arise in /3i^2i(t>f^ ,4'2,D, B) only from the corresponding 

terms c^_^ i('?^'i^^)^^~^V2 and c^_^ ^{<j}-^^)^^~^^ <j)"2^ . Assertion (2) now follows from 
the recursion relation 

— -^c'^ and — ^ ^± 

f-^,! — ^+l''^_2,l "^"^ ^^,2 ~ ^+l^£-2,2 • 

To prove Assertion (3), we first take Dirichlet boundary conditions. Let ^ > 4. 

Let <p\ \dM = for 7^ 1. No information is garnered concerning or . The 
term e7o arises in /3£_2((^^^'' + E(l)i,(j)2, D,B) only from the monomial 

+ -E'0i)(^"^V2- It now follows that 

e^"2 = (^-2)||Tc7-2,i = (^-2)2^- 
Since the coefficient c7_2 2 = 0' the term (j)''^^ (p^^ E^^~'^^ does not arise in the 
invariant;^/3£_2(</>i^^ + E(j)i,(j)2,D,Bg) and thus 
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Next we examine Neumann boundary conditions. We take 5 = and suppose 
<P^^\dM = for A; > 1. No information is garnered concerning e^g. Since c^^ = 0, 
the term e^^0i(?i2-E^^~^' and the term e^2^i'?^2^^^^^~^^ can arise in the invariant 
l3e-2{(i)x^+E(i)i,(i)2, D, B) only from the term 0^2(^1^^ +E(i)i)'^^-^'i We conclude 

The argument that = is similar and is therefore omitted. This establishes 

Assertion (3). 

To examine Assertion (4), we assume \qm = 'l>'\' \dM = 0- Again, we set £^ = 0. 
We study the terms ■^S(j)^-l^'^^ <j)2 and d^2'5'0i^~^V2^^- The case £ = 4 is a bit 
exceptional as these terms arise in h{4>i\4>2,D,B) only from 277" 1 5(</.fV 
and from 27r-i/2|S'(?!>f ')^i!>^^\ This shows that 

For £ > 6, these terms decouple and the recursion relation proceeds without com- 
plication to show 

~ ITT'^/-2,1 = ^"^^ — ITT'^^-2,2 ~ ■ '-' 

We can relate Neumann and Dirichlet boundary conditions. Let M := [0, 1] and 
let b e C°°{M). Let edr be the inward unit normal; £(0) = 1 and £(1) = — 1. 
Define: 

A:=dr + h, A*:=-dr + b, Di:=A*A, := AA* , 

S:=eb, B+:=eA, Ei:=b'-b^, E2 -.^ ~b' - b'^ . ^ 

Then = simply means A(I)\qm = 0. Furthermore Ei is the endomorphism 
defined by Di. 

Lemma 5.3. Adopt the notation established above. Let l>6 be even. 

(1) l3f^{^,,^2,D^,B+) = -j^Pe-2{MuM2,D2,B-). 

(2) e^-i =£-Ei, e+3 = (2 - e)Ee, d+^ = -2 • S^. 

Proof. Again, wc follow [12] to prove the first Assertion. Let {A^, 0^} be a complete 
spectral resolution of (_Di)g+. We obtain as above that 

We restrict henceforth to > since the contribution of zero eigenvalues to the 
above sum is zero. Let 



A, 

Then {A^, is a spectral resolution of D2 on Rangc(^) = ker(£)2)"'" with Dirichlet 
boundary conditions. Since A^j^laM = 0, the boundary terms vanish and we may 
express: 



'^'{Af)= [ (^/,V^)dvoU = ^ / (A/,A.^^)dvol„ 
= [ if, A*A,j>^) dvoU = ^A^T^ ^ (/) • 

H J M 
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This then permits us to express 

which yields the identity 

-dtPi^i ,ct>2,D^,B+) (i) = p{A<l>^ , Act>2 ,D2,B-) (t) . 

Assertion (1) now follows by equating terms in the asymptotic expansion in exactly 
the same fashion as was used to establish Assertion (1) of Lemma 5.2 (the extra 
negative sign can not be absorbed into D). 

We apply the relations of Equation (5.c) and use the fact that e^^ = cj_2 2 = 
to examine 

{<i>^^U2b^'-'\ <i>^^U2bb^'-'\ <t>^^u^^h^'-'\ <t>^^u^^\b^r-'^} . 

The assumption that £ > 6 is employed to ensure that S'^{(j){ '^^2^' + 4>i^4>2 ^'') 
does not produce such a term. We compute at the boundary component x = Q: 

el,<l>['U2Et''> = -Ee4'U2b('-'^ + 2 • E,4'U2bb('-^'^ + 

etA'^c^^^^Ef^ = elJ^U^M-') - elA'U^2\bT-'^ + 

d+,Sc^\'U2Et'^ = diAi^c^2bb^'-^) + 

-lh''7-2,AA^ + fo-^i (4'^ + bcf>2) + (4'^ + b^iM^ + 602)(^-^)} 

= -S,0f)^2&(^-^)-S,(£-2),/.(^V266(^-3)-2(£-2)S,<^W4'^^'('-') + .-., 

= -I^'=7-2,l{-</>l'V2&M^-3) _ ^(1)41)^,(^-3) _ ^(l)4l)(;,2)(.-4)} + .... 

This gives us the following relations: 

(a) <^W<^26(^-^) : -S, = -S,, 

{b) 4'^^'^cf>2bb^'-^) : 2 • + 43 = - 2) + j|Te7_2_i, 

(c) ^'^^U'i^b^'-^) : e+3 = -2{i - 2)S, + ^^ej_^^^, 
id) ##(6^)(^-4): -e+3 = ^e7_,,, . 

This then yields the following 3 relations: 

(1) (c)+(d): = -2(£-2)S,+2.j^e7_2,i so e7_2,i = (^-2)^-S, = (£-2)S,_2. 

(2) (d)-(c): -2e+3 = 2{i - 2)S^ so e+3 = (2 - 

(3) (c)-(b): -43 + e+3-2-S^ = -(£-2)S^ so ^+3 = e+3 + (^-4)S^ = -2•S^ □ 
We now work in dimension m > 2 to examine 

= / {ciAiU2 + e7i<Ai<^2^(^-'' + r7i<^i02P^-'^ + ...| dvol„_i . 

JdM ' 

Let Ml := [0, 1] and a € C°°{Mi) satisfy a\oMi = 0. Let 

Di := -d^, M2 := Ml x S\ D2 := Di - e'^^^''^^^ . 

Lemma 5.4. 

(1) Ife > 2, then = /3f^(l,e«M,-a2,S-). 

(2) r,- = i(£-2)S,. 



16 



M. VAN DEN BERG, PETER GILKEY, AND K. KIRSTEN 



Proof. We follow the treatment in [12] to prove Assertion (1). We consider the 
function u{r,t) = e~*^i'^~ 1. This solves the equations 

{dt+Di)u = 0, liinu{-,t) = 1 in L^{Mi), B-u = 0. 
Since u also solves the equations 

{dt+D2)u = 0, Yimu{-,t) = 1 in L^{M2), B-u = 0, 

we also have that u{-,t) = e~*^^'^l as well. Since cIvoIms = e"drd9, 

PM,{l,e-",D2,B-){t) = u{r,t)e-"^'-'>e"^'-'>dedr 

Jr=Q J 0=0 

= 27r / u(r,t)rfr = 27r^Mi(l,l,-Di,B-)(i). 

Jr=0 

Since the structures are flat on Mi, /3f (1, 1, Di,B-) = for ^ > and A|f^ 1=0. 
We equate terms in the asymptotic expansion to see jSf^'^ {1, e~°'^^\ D2, B~) = 
for £ > as well. 

We apply Assertion (1). We use the formalism of Equation (l.e). We have 
ds%i^ = dr^ + e^^^^^dO'^ where a(0) = and a{r) = near a = 1. We compute: 

ri22 = r2i2 = -r22i = e^"a^^\ wi = ie-2«r22i = -^a(^\ 
002 = 0, £(^-2) = iaW + 

^^'=0 + ..., <^W = -aW + ..., 

We examine the coefficient of a^^^ in j3i for i even: 

el,<Pi<j>2E('-') = ^e-EeaW + ...., 

r70i02Pmm^^ = -r^aW + ... . 

It now follows from Assertion (1) that rj = ^{i — 2)S^. This completes the proof 
of Lemma 5.4 and thereby completes the proof of Theorem 1.9 as well. □ 



6. Estimating the heat trace asymptotics on a closed manifold 

In this section, we shall prove Theorem 1.2. We shall proceed purely formally 
and shall use the discussion in Sections 1.7-1.8 of [27] (which is based on the Seeley 
calculus [44, 45]) to justify our formal procedures. As in Equation (l.b), let 

D = -g'^d^.d^^ - A^'d^^ - B 

be an operator of Laplace type. Throughout this section, C = C{M, g, D) will de- 
note a generic constant which depends only on (M, g, D) (and hence also implicitly 
on m) but not on n; c(m) will denote a generic constant which only depends on m. 
If we take D = A^, then C = C{M, g). 

We introduce coordinates ^ = (^1, ■■■,£,m) on the cotangent bundle to express a 
covector in the form ^ = ^jda;'. The symbol of D is P2{x, ^) +pi{x, ^) +Po{x) where: 

P2ix, := g'^{x)£,i^,j, pi{x, := A''{x)£,k, and pa = B . 

There are suitable normalizing constants involving factors of \/— T which we ignore 
in the interests of simplicity henceforth since they play no role in the estimates we 
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shall be deriving. Let C := C — [0,00) be the slit complex plane and let A € C. 
Following the discussion in Lemma 1.7.2 of [27], one defines inductively: 

rn{x,^,X)--=-ro{x,^,X)- ^ d^pkix,^) ■ d^rjix,^,X)/a\ . (6.a) 

\a\-\-j+2—-k—n.j<n 

In this sum fc = 0, 1, 2 and |q;| < 2 — k. The symbol of e~*^ is given by: 

eo{x, (,,t) + ... + e„(x, ^, t) + ... 
where, following Equation (1.8.4) of [27], one sets: 

en{x,^,t) := / e"*-^r„(a;,^, A)dA; 

iTTV— 1 J J 

here 7 is a suitable contour about the positive real axis in the complex plane. Then, 
following Equation (1.8.3) of [27], one may obtain the local heat trace invariants of 
Equation (l.c) by setting: 

-1 



anix,D) = \^^^det{g^j)j J e„(x,^,l)d^. (6.b) 

To measure the degree of an expression in the derivatives of the symbol, we set: 

degree(d^(?*^') = |a|, degree(d;;!A'=) ^ \a\ + 1, dcgree(d;;!B) = |a| + 2 . 

Note that if D is the scalar Laplacian, then B = and A'^ = g~^dxig^-'g has degree 
1 in the derivatives of the metric so this present definition is consistent with our 

previous definition in this special case. It is immediate from the definition that ro 
is of total degree in the jets of the symbol of D. Furthermore, since 

degrec((i|pfe) = 2 — A: and degvee{d'^r j) = \a\ + degree(rj) , 

we have by induction that 

degree(r„) = n . (6.c) 
There is a similar grading on the variables (^, A). One defines: 

weight(^i) = 1 and weight(A) = 2 . 

It is then immediate that rg has weight —2 in (^, A). Clearly 

weight (d^^pfe) = fc — |a| and weight(d^rj) = weight (rj) . 

Thus it then also follows by induction from Equation (6. a) that 

weight(r„) = — 2 — n . (6.d) 

Let n be odd. Since the weight of r„(x,^, A) is —n — 2 in (^, A), it follows that 
e-n{x, 1) is an odd function of ^ and hence the integral in Equation (6.b) vanishes 
in this instance. This yields an{x,D) = for n odd. Let [•] be the greatest integer 
function. 

Lemma 6.1. 

(1) We may expand Vn in the form: 

2n+l 

rn{x,^,X)= ^ qn,m,j,i3{x,g)^'^rj,{x,^,X) . 

i=[in] + l |/i|=2j-«-2 

(2) There exists a constant C{M,g) so that ifn = 2n > and if\P\ = 2j — n—2, 
then 

C{M,gr 



e-^rl{x,$„X)fdXd£, 



< 
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Proof. We apply the recursive scheme of Equation (6. a) to obtain an expression for 
r„ of the form given in Assertion (1). By Equation (6.c), r„ has degree n in the 
derivatives of the symbol of D. Thus there arc at most n x-derivativcs of ro which 
are involved in the process. Each x-derivative of ro adds one power of ro (other 
variables can be differentiated as well of course so we are obtaining an upper bound 
not a sharp estimate). Each step in the induction process adds 1 power of ro. Thus 
j <2n+ 1. By Equation (6.d), r„ is homogeneous of weight — n — 2 in A). Since 
\(3\ — 2j = —n — 2 and |/3| > 0, we may conclude that j > 1 + |n > [^n] + 1. 
Assertion (1) now follows. 

We use the Cauchy integral formula to estimate: 

1 



< 



/ 



ed^ 



The quadratic form g^^ is positive definite. Thus we may estimate > e\^\l for 
some £ = e{M,g) > where ~ £,i + ... + is the usual Euclidean length. 

Note that \^^\ < Since e^'^'^ < e^^l^l=, we may use spherical coordinates to 

estimate: 

C30 



7 



< 



1 



poo 

/ e-^'-Vl^l+'"drvol„_i(S'"-\3s— • 

Jr=0 



Since < 2j < 4n + 4 is uniformly and linearly bounded in n, we may rescale to 
remove s in e~^'" at the cost of introducing a suitable multiplicative constant. We 
may then evaluate the integral to estimate: 

- xr^^dxdi 



7 



<C{M,gy 



2 /• 



(.?-!)! 



Since j — 1 — ^\(3\ = n the desired estimate follows; the shift by m can be absorbed 
into C{M,g)'^ since we have restricted to n > 0. □ 

Let Df c C™ be the complex polydisk of radius e of real dimension 2m about 
the origin in C™ given by setting: 

D'^ ■.= {z = {zi,...,Zm)&C'" ■\zi\<e for 1 < i < m} . 

We let Df = of n R"* be the corresponding real polydisk. We also consider the 
submanifold of real dimension m in C™ (which is not the boundary either of the 
complex polydisk Df or of the real polydisk Df) given by: 

5e := C™ : \zi\ =e for 1 < i < m} . 

We consider the holomorphic m-form 

dw = {2'K^/^-V)~™'d'W\...dwm ■ 

Let / be a holomorphic function on the interior of which extends continuously 
to aU of D'f and let a is a multi-index. If z belongs to the interior of the polydisk 
DS:, then we shall define: 



f{w){wi - Zi} 



-1 — ai 



• (^m ) 



"■dw . 



wes. 



We may then use the Cauchy integral formula to represent: 
d^f{z) = a\lM) for zeint(DC). 
Let /? = P{i, a) be the multi-index {ai, .... ai^i,ai + 1, a^+i, . 

We introduce variables {fv} for the {g^^,A'',B} variables; we have a total of 
^m{m — 1) + m + 1 such variables. Since we are in the real analytic setting, we can 



, ctm)- We then have: 
(6.e) 
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choose real analytic coordinates about each point P of M which are real analytically 
equivalent to the polydisk D^{P) of radius 2 in such a way that the variables {f,^} 
extend continuously to D2{P) with fi, holomorphic on the interior of -D2 The 
functions \fy\ are uniformly bounded on D^{P). If G (P) and |w| € S^{P), 
then \zi — Wi\ > 1 and thus we have uniform estimates 

\Ia{U){z)\<C{M,D) for any u,a. (6.f) 

We decompose r„ in terms of monomials of the form 

rif ■ g'^'^ ■ ... ■ 9'^'^ ■ laAM ■ ... ■ . (6.g) 

Here wc assume degree{9^/^^} > since we have made explicit the dependence on 
the variables of degree 0. Thus b < n since, by Equation (6.c), r„ is homogeneous 
of degree n in the jets of the symbol. There are no 5'^ variables in ro- Each 
multiplication by d"p2 can add at most one g*^ variable; each multiplication by d^.pi 
or pq adds no g^^ variable. Each application of to rj does not add a g^^ variable 
(and can in fact reduce the number of g'^^ variables if they are differentiated). Thus 
the number of g^^ variables is at most n. Thus in considering monomials of the form 
given in Equation (6.g), we may assume a < n. Wc summarize these constraints: 

j<2n+l, -n- 2 = - 2j, a < n, and b<n. (6.h) 

Lemma 6.2. Let c{m) := 50m^. We can decompose r„ as the sum of at most 
c(m)"'n! monomials of the form given in Equation (6.g) satisfying the constraints 
0/ Equation (6.h) where the coefficient of each monomial has absolute value at most 
1. 

Proof. Since ro can be written as a single monomial with coefficient 1, we proceed 
by induction. 

(1) Consider —r^d^^p^ ■ dx^fn-i. Each k generates m terms so there are 
terms generated in this way. Differentiating Tq generates at most 3n terms 
since j < 3n by Equation (6.h). Differentiating the f/' variables generates 
at most n terms since a < n. Differentiating the I variables generates at 
most 6 + X] < 2n terms by Equation (6.e). Thus we generate at most 
m^(3n + n + 2n) — 6m^n terms from each monomial of r„_i. This can be 
written in terms of at most 

6m^n • c(m)™^^(n — 1)! = 6TO^c(m)™^^n! monomials. 

(2) Consider — ro9{j.^ ^^k^P"^ ' ^^k^ ^xk^ rn-2 ■ A similar argument shows this gen- 
erates at most m^(6n)(6(n — 1)) new terms from each monomial of rn-2. 

This can be written in terms of at most 

36TO^n(n — 1) • c(m)"~^(n — 2)! < 36m^ • c{m)"~^nl monomials. 

(3) Consider —rfjA^^^rn^i. This can be written in terms of at most 

m ■ c{m)"^~^{n — 1)! < m^c{m)"^~^nl monomials. 

(4) Consider —roA'^dxf,rn-2- This can be written in terms of at most 

6mn . c(m)"^^(n — 2)! < 6m^c(m)"^^n! monomials. 

(5) Consider —roBrn-2- This can be written in terms of at most 

c(m)""^(n - 2)! < m'^c{m)'^~^n\ terms. 

The above argument shows that r„ can be decomposed as the sum of at most of 
50to^ • c(m)"~^n! = c(m)" • n\ monomials each of which has a coefficient of absolute 
value at most 1. □ 
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Proof of Theorem 1.2 (1). We consider monomials where the coefficient has abso- 
lute value at most 1. We have shown that there exists a constant c(m) so that r„ 
can be written in terms of at most c(m)"n! such monomials. We may then use 
the constraints of Equation (6.h), the estimates of Equation (6.f), and the esti- 
mate of Lemma 6.1 to construct a new constant C{M,g) and complete the proof 
of Theorem 1.2 (1) by bounding: 

\an{x,D)\ < c(m)"n! • C{M,g,Df^ ■ C{M,gr^ < C{M,g,DYn\. □ 

Proof of Theorem 1.2 (2). Let P be a point of a closed real analytic Riemannian 

manifold {M,g). Let / be a real analytic function on M so that df{P) ^ 0. Since 
/ is continuous and M is compact, |/| is bounded. By rescaling and shifting /, 
we may suppose without loss of generality that /(P) = and that |/(a;)| < 1 for 
all points x of M. We make a real analytic change of coordinates to assume that 
g^^{P) = 6ij and that f{x) = Cf ■ xi near P. We shall choose = ±1 recursively 
and define: 

oo 

h{x)=J2^k2-''f{xf'. 

k=3 

This series converges uniformly in the real analytic topology so h is real analytic. 
Let £n{-) be a generic invariant which only depends on the parameters indicated. 
Let gh = e^'\g. Let n > 3. We use Theorem 1.8 to see that: 

(d^^ym = e^2-«cf (2n)! +f i(ci. 

Tg^{P) ~ Cm{d'^^h){P) + lower order terms for some \cm\ > 1, 

a2n{P, Ag) = (-1)"-! (2|fl)T^""^^ + lower order terms 
= Cmj^^cfen2-''{2ny.+£i{ei,...,en-i,g). 

We set 

_ f -1-1 if Cm£U^l,:.,£n-l,9)>0] 
\ -1 if Cm£U^l,-,£n-l,9) <0 J 

With this choice of en, there is no cancellation. As ^ ^-"^^ > ^ for n > 3, we 
obtain the desired estimate: 

|a2^(P,A,)| >c„(Jfi^cf2-'^(2n)!> ^cf2-«.n!> (Ac2)"n!. □ 

7. Growth of heat content asymptotics 

This section is devoted to the proof of Theorem 1.4. We first examine a product 
manifold [0, 1] x N. Let {e^} be a sequence of signs to be chosen recursively. We 
replace the function f{x) of the previous section by sin(x) and define: 

oo 

h{x) ■.= ^e^2-''sm{xf'' . 
i/=i 

This scries converges in the real analytic topology to a real analytic function h 
which is periodic with period 27r and which satisfies h{0) = /i(27r) = 0. We set 

gM := e^^'{dx^ +gN) ■ 

The inward unit normal is given at by r/(0) = dx and at 27r by i/(27r) = —dx. If 
j is odd, then {dih}{0) = {9^/i}(27r) = since h is an even function. And clearly 
we have that {{dl)h}{0) = {{—dx)-'h}{2'jT) if j is even. Consequently 

/iW(0) = /i(-'''(27r) for any j. 
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This ensures that the behaviour of h is the same on the boundary components and 
gives rise to the factor of 2 vo1to-i(-^'^, ffiv) in Equation (7.a) below. We have: 

Since rn > 2, there is a non-zero constant Cm with \cm\ > 1 which only depends on 
m and not on i so that: 

Wc may then apply Theorem 1.9 to express: 

/3|f (1, 1, AM,g„,S-) = {i(2Z- 2)S2^-c„2-^"(2£)! • 2 voU_i(iV,3jv)} 

■'2 A 



•••,£^_l,fi'iv) • 



(7.a) 



Set 



+1 if £|.-(ei, ....e^.ijgjv) > 



2e 

-1 if f|^-(£i,...,e^-_i,5jv) < J 
Since there is no cancellation in Equation (7. a), we may estimate: 

l^ff (1, 1, Am,sm,S-)| > ^{2£- 2)'E^iCm{2l)\ei2-~^ ■ 2vol„_i(Ar,5w) 
The desired estimate in Assertion (1) of Theorem 1.4 now follows since: 



^(2^-2)S2,-c„(2/)!£,-2- 



> (21 - 2)^—...--^2-h ■ 2 ■ 3... • 2£ 



= f^^2 • 4 • ... • 2^ > —2^e\ > l\ for £ > 3. 
2^+1 - 14 - 

We now turn to the case of the ball and apply a similar analysis to establish 
Assertion (2) of Theorem 1.4. The functions sin(a;) is now replaced by the function 
{x\ + ... + a;^ — 1)^'^, the operator dx is replaced by the radial derivative 9^, and 
the boundary components a; = and x = 2-k are replaced by the single boundary 
component r = 1. The remainder of the argument is the same and is therefore 
omitted. □ 



Acknowledgments 

Research of M. van den Berg partially supported by the London Mathemat- 
ical Society Grant 41028. Research of P. Gilkey partially supported by Project 
MTM2009-07756 (Spain), by DFG PI 158/4-6 (Germany), and by Project 174012 
(Serbia). Research of K. Kirsten supported by National Science Foundation Grant 
PHY-0757791. It is a pleasure to acknowledge useful conversations with Professor 
Michael Berry (University of Bristol) concerning these matters. 

References 

[1] p. Amsterdamski, A. Berkin, and D. O'Connor, "bg 'Hamidew' coefficient for a scalar field". 

Classical Quantum Gravity 6 (1989), 1981-1991. 
[2] I. Avramidi, "Covariant methods of studying the nonlocal structure of an effective action" , 

Sov. J. Nud. Phys. 49 (1989), 735-739. 
[3] 1. Avramidi, "The nonlocal structure of the one-loop effective action via partial summation 

of the asymptotic expansion", Physics Letters B 236 (1990), 443-449. 
[4] I. Avramidi, "A covariant technique for the calculation of the one-loop effective action" , Nucl. 

Phys. B 355 (1991), 712-754. 
[5] M. van den Berg, "The Casimir effect in two dimensions". Physics Letters A 81 (1980), 

88-90. 

[6] M. van den Berg, "On finite volume corrections to the equation of state of a free Bose gas" , 

Helvetica Physica Acta 56 (1983), 1151-1157. 
[7] M. van den Berg, "Asymptotics of the heat exchange", J. Fund. Anal. 206 (2004), 379-390. 



22 M. VAN DEN BERG, PETER GILKEY, AND K. KIRSTEN 

[8] M. van den Berg, "On the free boson gas in a weak external potential", Physics Letters A 
78 (1981), 219-222. 

[9] M. van den Berg, and E. B. Davics, "Heat flow out of regions in K""', Math. Z. 202 (1989), 
463-482. 

[10] M. van den Berg, and J.-F. Le Gall, "Mean curvature and the heat equation", Math. Z. 215 
(1994), 437-464. 

[11] M. van den Berg, and S. Srisatkunarajah, "Heat flow and Brownian motion for a region in 

]R2 with a polygonal boundary", Probab. Th. Rel. Fields 86 (1990), 41-52. 
[12] M. van den Berg, S. Desjardins, and P. Gilkey, "Functorality and heat content asymptotics 

for operators of Laplace type", Topol. Methods Nonlinear Anal. 2 (1993), 147—162. 
[13] M. van den Berg, and P. Gilkey, "Heat content asymptotics of a Riemannian manifold with 

boundary", J. Funct. Anal. 120 (1994), 48-71. 
[14] M. van den Berg, P. Gilkey, A. Grigoryan, and K. Kirsten, "Hardy inequality and heat 

semigroup estimates for Riemannian manifolds with singular data" , 

http://arxiv.org/abs/1011.1726. 
[15] M. van den Berg, P. Gilkey, K. Kirsten, and V. A. Kozlov, "Heat content asymptotics for 

Riemannian manifolds with Zaremba boundary conditions", Potential Analysis 26 (2007), 

225-254. 

[16] M. van den Berg, P. Gilkey, and R. Secley, "Heat content asymptotics with singular initial 

temperature distributions", J. Funct. Anal. 254 (2008), 3093-3122. 
[17] M. V. Berry, and C. J. Howls, "High orders of the Weyl expansion for quantum billiards: 

resurgence of periodic orbits and the Stokes phenomenon", Proa. R. Sac. Land. A 447 (1994), 

527-555. 

[18] M. Bordag, G. L. Klimchitskaya, U. Mohideen, and V. M. Vostepanenko, "Advances in the 
Casimir effect", International Series of Monographs on Physics 145. Oxford University Press 
(2009). 

[19] T. Branson, P. Gilkey, K. Kirsten, and D. Vassilevicli, "Heat kernel asymptotics with mixed 

boundary conditions", Nucl. Phys. B 563 (1999), 603-626. 
[20] T. Branson, P. Gilkey, and B. 0rsted, "Leading terms in the heat invariants", Proc. Amer. 

Math. Soc. 109 (1990), 437-450. 
[21] R. Brooks, P. Perry, and P. Yang, "Isospectral sets of conformally equivalent metrics", Duke 

Math. J. 58 (1989), 131-150. 
[22] H. S. Carslaw, and J. C. Jaeger, Conduction of Heat in Solids, Clarendon Press, Oxford, 

(2000). 

[23] A. Chang, and P. Yang, "Compactness of isospectral conformal metrics on S^" Comment. 

Math. Helv. 64 (1989), 363-374. 
[24] L. Geisinger, and T. Weidl, "Universal bounds for traces of the Dirichlet Laplace operator" , 

J. Land. Math. Soc. 82 (2010), 395-419. 
[25] P. Gilkey, "Recursion relations and the asymptotic behavior of the eigenvalues of the Lapla- 

cian", Compositio Math. 38 (1979), 201-240. 
[26] P. Gilkey, "Leading terms in the asymptotics of the heat equation" , Geometry of Random 

Motion, Ithaca 1987, Contemp. Math. 73, 79-85, Amer. Math. Soc, Providence, R. L, (1988). 
[27] P. Gilkey, Invariance Theory, the heat equation, and the Atiyah-Singer index theorem, CRC 

Press Boca Raton, Fl., (1994). 
[28] P. Gilkey, Asymptotic Formulae in Spectral Geometry, Stud. Adv. Math., Chapman & 

Hall/CRC, Boca Raton, FL, (2004). 
[29] P. Grcincr, "An asymptotic expansion for the heat equation", 1970 Global Analysis Proc. 

Sympos. Pure Math. XVI Berkeley, Calif., 133-135, Amer. Math. Soc, Providence, R. I., 

(1968). 

[30] P. Greiner, "An asymptotic expansion for the heat equation" , Arch. Rational Mech. Anal. 
41 (1971), 163-218. 

[31] C. H. Howls, and S. A. Trasler, "High orders of Weyl series; Resurgence for odd balls", J. 

Phys. A: Math. Gen. 32 (1999), 1487-1506. 
[32] K. Kirsten, "The as heat kernel coefficient on a manifold with boundary", Classical Quantum 

Gravity, 15 (1998), L5-L12. 
[33] K. Kirsten, "Spectral Functions in Mathematical Physics", Chapman and Hall/CRC, Boca 

Raton, Fl., (2002). 

[34] T. Krainer, "On the expansion of the resolvent for elliptic boundary conta<;t problems", Ann. 

Global Anal. Geom. 35 (2009), 345-361. 
[35] D. M. McAvity, "Surface energy from heat content asymptotics", J. Phys. A: Math. Gen. 

26 (1993), 823-830. 

[36] P. McDonald, and R. Meyers, "Dirichlet spectrum and heat content", J. Funct. Anal. 200 
(2003), 150-159. 



GROWTH OF HEAT TRACE AND HEAT CONTENT ASYMPTOTIC COEFFICIENTS 23 



[37] H. P. McKcan, and I. M. Singer, "Curvature and the eigenvalues of the Laplaciaii'' , J. Dif- 
ferential Geom. 1 (1967), 43-69. 

[38] S. Minakshisundaram, and A. Pleijel, "Some properties of the eigenfunetions of the Laplace 
operator on Rieinannian manifolds", Canad. J. Math. 1 (1949), 242-256. 

[39] B. Osgood, R. Phillips, and P. Sarnak, "Compact isospectral sets of surfaces" , J. Fund. Anal. 
80 (1988), 212-234. 

[40] C. G. Phillips, and K. M. Jansons, "The short-time transient of diffusion outside a conducting 
body", Proc. R. Soc. Land. A 428 (1990), 431-449. 

[41] Yu. Safarov, and D. Vassilicv, The asymptotic distribution of eigenvalues of partial differ- 
ential operators. Translated from the Russian manuscript by the authors. Translations of 
Mathematical Monographs 155. American Mathematical Society, Providence, R. I., (1997) . 

[42] A. Savo, "Uniform estimates and the whole asymptotic series of the heat content on mani- 
folds", Geom. Dedicata 73 (1998), 181 214. 

[43] R. Seeley, "Singular integrals and boundary value problems", Amer. J. Math. 88 (1966), 
781-809. 

[44] R. Seeley, "Complex powers of an elliptic operator" , Singular Integrals, Proc. Sympos. Pure 
Math., Chicago, 111., 288-307. Amer. Math. Soc, Providence, R.I. (1966). 

[45] R. Seeley, "Topics in pseudo-differential operators", 1969 Pseudo-Diff. Operators (C.I.M.E., 
Stresa, 1968), 167-305. Edizioni Cremonese, Rome. See also Pitman Research Notes in Math- 
ematics Series, 359. Longman, Harlow, 1996. 

[46] R. Seeley, "The resolvent of an elliptic boundary value problem", Amer. J. Math. 91 (1969), 
889-920. 

[47] R. Seeley, "Analytic extension of the trace associated with elliptic boundary problems", 

Amer. J. Math. 91 (1969), 963-983. 
[48] I. Travenec, and L. Samaj, "High orders of Weyl series for the heat content", 

arXiv:1103.0158vl. 

MvdB: School of Mathematics, University of Bristol, University Walk, Bristol BS8 
ITW, UK 

E-MAIL: M.VANDEXBERKCQ'BRIS.AC.LIK 

PC: Mathematics Department, University of Oregon, Eugene OR 97403 
E-mail: gilkey@uoregon.edu 



KK: Department of Mathematics, Baylor University, Waco, Texas, TX 76798, USA 
E-mail: Klaus_Kirsten@baylor.edu 



